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ABSTRACT 

We make a simple analytical study of radial profiles of dark matter structures, with 
special attention to the question of the central radial density profile. We let our the- 
oretical assumptions be guided by results from numerical simulations, and show that 
at any radius where both the radial density profile, p, and the phase-space-like density 
profile, p/cr £ , are exact power laws, the only allowed density slopes in agreement with 
the spherical symmetric and isotropic Jeans equation are in the range 1 < j3 < 3, 
where [3 = —d\np/d\nr. We also allow for a radial variation of these power laws, as 
well as anisotropy, and show how this allows for more shallow central slopes. 
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1 INTRODUCTION 

The formation and evolution of Dark Matter (DM) struc- 
tures is in principle very simple since it only involves gravity. 
Despite this fact, density profiles of dark matter halos have 
become one of the most challenging issues for our under- 
standing of cold dark matter structure formation. Numeri- 
cal simulations provide predictions of steep central density 
cusps with power law slopes, p ~ r -13 , with f3 from 1 to f .5 
within a few percent of the virial radius of the halo (Navarro, 
Frenk & White 1996; Moore et al. 1998). Recent careful 
studies (Diemand, Moore & Stadel 2004; Reed et al. 2003; 
Navarro et al. 2004) indicate that the resolved region has 
still not converged on a central density slope, so in principle 
the central power slope may be even shallower. 

The steep inner numerically resolved slopes are, how- 
ever, not supported by observations. By measuring the ro- 
tation curve of a galaxy one can in principle determine the 
density profile of its DM halo. Low surface brightness galax- 
ies and spirals, where the observed dynamics should be DM 
dominated, seem to show slowly rising rotation curves (Ru- 
bin et al. 1985; Courteau 1997; Palunas & Williams 2000; de 
Blok et al. 2001; de Blok, Bosma & McGaugh 2003; Salucci 
2001; Swaters et al. 2002; Corbelli 2003) indicating that 
these DM halos have constant density cores. Galaxy clusters, 
where baryons can play even less of a role, may show a sim- 
ilar discrepancy. Arcs (Sand, Treu & Ellis 2002) and strong 
lensing fits of multiple image configurations and brightnesses 
(Tyson, Kochanski & dell'Antonio 1998) also indicate shal- 
low cores in clusters. All these observations could be in 
agreement with N-body simulations only if either the very 
central region is really not cuspy, or if cusps could somehow 
be erased during galaxy formation. 

It is therefore very important to understand if the pure 



dark matter central density slopes can really be as steep as 
indicated by numerical simulations, in order to understand 
if one needs to invoke baryonic physics to reach agreement 
with observations. Baryonic structures are often observed 
to to have central cores, which is possibly even understood 
theoretically (Hansen & Stadel 2004). 

Several attempts have been made for an analytical 
derivation of the density profile (Bertschinger 1985; Syer 
& White 1998; Subramanian, Cen & Ostriker 2000; Hiotelis 
2002; Dekel et al. 2003; Manrique et al. 2002), and none seem 
to present a clear and simple explanation for the findings of 
N-body codes. We will not attempt to answer the very diffi- 
cult question of the actual formation of DM structures here, 
but will instead simply consider the Jeans equation and ask 
which equilibrium DM structures are in agreement with the 
spherically symmetric Jeans equation. We will be guided by 
the findings of numerical simulations and only consider the 
special cases where the phase-space-like density, p/cr E , is a 
power law in radius for some positive e. 

The normal use of the Jeans equation for collisionless 
systems (Hernquist 1990; Tremaine et al. 1994) is to as- 
sume a given radial density profile, p(r), and then solve the 
Jeans equation to get the corresponding velocity dispersion, 
<r 2 (r). This can be done analytically for sufficiently nice den- 
sity profiles, and can always be done numerically. The ba- 
sic result is that the Jeans equation can allow for almost 
any shape of the density profile. An alternative approach 
is instead to assume the form of the phase-space density, 
p/a 3 (r), and then solve the Jeans equation to get the corre- 
sponding density profile (Taylor & Navarro 2001). Also this 
can always be done numerically, and even analytically for 
sufficiently nicely behaving velocity dispersions. 

We will show below that for sufficiently simple phase- 
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space (like) densities, this approach can provide analytical 
insight into the allowed range of density profiles. One ex- 
ample hereof is that if both the central density profile and 
the phase-space-like density are exact power laws, then the 
central density profile of an isotropic DM structure cannot 
be more shallow than an NFW profile with (3 = 1. 



2 EXACT POWER LAWS 

Let us first consider the case where the coarse grained radial 
density profile is an exact power law 



p ~ r 



(1) 



at a given radius. Now, Taylor & Navarro (200f ) observed 
that the phase-space density, p/a" 3 , from N-body simulations 
approximately follows a power law. Recent high resolution 
N-body results confirm that this is approximately correct in 
the equilibrated inner region, where substructures are unim- 
portant (Diemand, private communication). * We will here 
make a slightly weaker assumption, namely that a phase- 
space- like density profile is an exact power law in the very 
central region 

£ ~ r"« , (2) 

with unknown real numbers e > and a. Taylor & 
Navarro (2001) found e — 3 and a — 1.875. Recent high res- 
olution N-body simulations do infact support this assump- 
tion for the very central numerically resolved region with e 
of the order 2-3. We will here not attempt to understand 
why the phase-space-like density is a power law in the cen- 
trally resolved region. A reason therefore must be sought at 
a deeper level, maybe through a solution to the collisionless 
Boltzmann equation. 

For spherically symmetric and isotropic systems the 
Jeans equation can be written (Binney & Tremaine 1987; 
Taylor & Navarro 2001) through the use of Poisson equa- 
tion (for a self gravitating system) 



d_ 

dr 



-r 2 d 
Gp dr 



(c 2 p) ) = 47rpr 2 



(3) 



This assumption of a spherical, isotropic system is guided by 
the numerical N-body results in the central part of the DM 
structure. The inclusion of anisotropy, e.g. Ap = 1 — Vg/v%, 
gives another term in the Jeans equation, 2Apd(ra 2 )/dr, 

and will therefore increase the space of solutions < , and we 
will later show how. We will leave non-spherical structures 
for a later analysis. 

Under the assumption of power laws, eq. (3) can now 
be written 



-CiC 2 r 2(a -« /e = C 3 r a - 



(4) 



* Taylor & Navarro (2001) considered spherical bin averages of 
p(r) and c 3 (r), and then took the ratio, p/cr 3 . The actual phase- 
space density, which is spherical averages of p/cr 3 , differs due to 
substructures (Diemand, Moore & Stadcl 2004; Arad, Dekel & 
Klypin 2004). We consider only the equilibrated inner region of 
the DM structure where there is no difference, 
t We use an unusual notation for the anisotropy parameter, Ap, 
to avoid confusion with the (3 in the density profile. 



where the two coefficients Ci and Ci come from the radial 
differentiations, e.g. Ci = dln(<r 2 p)/cilnr, and the last coeffi- 
cient Cz is a positive constant. Clearly, the radial power- laws 
in eq. (4) have to agree, giving 



(3 



2(e-a) 



(5) 



(e-2) • 

Moreover, for the Jeans equation, eq. (3), to make sense, 
the product C\ C2 must be negative. If the product, Ci C2, 
is positive, then there will be something negative on the lhs 
of eq. (4) and something positive on the rhs. We can thus 
find the points where the Jeans equation breaks down by 
solving Ci C2 = 0. This is a simple quadratic equation in a, 
with solutions 



a = 2 ± (e/2 - 1) 



(6) 



That is, when a has the value in eq. (6) then the lhs of the 
Jeans equation is zero. 

Combining the two results in eqs. (5, 6) tell us, that the 
only allowed values for the density slope are in the range 



1 < P < 3, 



(7) 



and one thus concludes that for this most simple case of pure 
power laws, the central density profile of pure dark matter 
structures cannot be more shallow than p ~ r _1 . Please note 
that this result is obtained for rather general power laws like 
eq. (2) with any value of e and a. The results of recent N- 
body simulations tell us that locally the density profile can 
be approximated by a power law, and furthermore one can 
always find a value for e such that eq. (2) holds true locally. 
Therefore, for any radius in the resolved region of N-body 
simulations where substructure is not important, the density 
profile must be in the range of eq. (7) . 

When one includes non-isotropic systems, where Ap 7^ 
0, then one finds the lower limit to be 



Anin = 1 + Ap , 



(8) 



which implies that for a negative Ap one can have more 
shallow profiles, e.g. a core in density for sufficiently circular 
orbits. For purely radial motion the most shallow profile is 
/3 m in = 2, which is in agreement with numerical findings 
for the spherical infall model (Lokas & Hoffman 2000). In 
general Ap can take any value in the range, —00 < Ap < 1 
(Mamon & Lokas 2004), naturally constrained by f3 being 
always non-negative (Tremaine et al. 1994; Hansen et al. 
2004), however, it should be kept in mind that numerical 
simulations find very little anisotropy in the very central 
region, Ap ~ (Moore et al. 2001). The central isotropy can 
also be understood from a fundamental statistical mechanics 
point of view (Hansen et al. 2004) . 



3 ALMOST POWER LAWS 

If the density profile and the phase-space-like density profile 
are not exact power laws, then our findings in eq. (7) may 
potentially be invalid. In order to investigate this question 
we can make expansions around power laws. As we will see, 
this will also indicate how far beyond the resolved region it 
makes sense to extend our findings. We therefore write 



-P(r) 



(9) 
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(10) 



£_ ^ .-«(>■) 

(7« 

where /3(r) and a(r) are now slowly varying functions of 
radius. We choose a fixed, radially-independent e. Now, let 
us make a Taylor expansion around the radius r_i, where 
[3=1, using (3' = d(3/d\nr and a' = da/dinr. It should 
therefore be kept in mind that this Taylor expansion only 
holds sufficiently nearby the point of expansion, such that 
the higher derivative can be ignored, (3' S> (3" lnr. The Jeans 
equation again looks like eq. (4) , and the radial powers again 
lead to the expression for (3 in eq. (5). However, the coeffi- 
cients are now different 

Ci = -13 + 2 {a - P) + ln(r/r-i) - (a' - (3' (1 + e/2)) 



C 2 = 



2 (a - /3) 2 , , rflnd 

1 H h ln(r/r_i) - (a - /? J + -= 

e e v y dlnr 



where ln(r/r_i) appears since we make the expansion 
around r_i. The first coefficient, Ci, is the one determin- 
ing the most shallow slope, and we find that one has Ci = 
when 

a = | + 1 - ^ ln(r/r_i) (a' - + e/2)) (11) 
which through eq. (5) implies 

Anin = 1 + - ln(r/r-i) (a' - /?' (1 + e/2)) . (12) 

One sees that the density slope within radius r_i can be 
slightly more shallow than f3 = 1. This is in agreement with 
the numerical findings of Taylor & Navarro (2001). 

The most recent simulations indeed seem to indicate 
that the phase-space-like density profile, eq. (10), can indeed 
be well fitted with a power law in the central resolved region, 
where e is found to be in the range e ~ 2 — 3 (Diemand, pri- 
vate communication). One can therefore always choose the 
epsilon in such a way that a' = locally. Furthermore, it is 
interesting to make a comparison with a recent beautiful fit- 
ting formula valid for the entire resolved region, as presented 
in Navarro et al. (2004) 



0n(t) = - 



dlnp 
dlnr 



r 

r-i 



(13) 



where r_i is the radius where (3 = 1. Since this formula gives 
(3 smaller than 1 it may lead to a constraint on e from the 
phase-space- like density. If this formula, eq. (13), is consis- 
tent with the spherical and isotropic Jeans equation, then 
this /3n must be larger than the smallest allowed (3 as de- 
termined from eq. (12), in the range where the Taylor ex- 
pansion leading to eq. (12) is valid, i.e. in the vicinity inside 
r_i. This is solved by 



e > 2. 



(14) 



When the numerical N-body simulations reach the level 
of resolution where they can resolve inside r_i, it will be 
straightforward to test the validity of the Jeans equation 
through the phase-space-like density profile. Thus, if one 
finds numerically that the phase-space-like density is indeed 
a power law, and only so with e < 2, while simultaneously 
the density profile is sufficiently close to a power law (to 
assure validity of the Taylor expansion, as quantified after 
eq. (10)), then either that region is not resolved numerically, 
or the formula (13) breaks down. Clearly, it is possible that 



the simulations will find no density slope more shallow than 
(3 = 1, in which case there is no constraint on e. 



4 CONCLUSIONS 

Recent numerical dark matter simulations show that the 
phase-space- like density profile, p/o 1 , is well fitted locally 
with a simple power law with e of the order 2-3. We show 
that when the radial density profile is an exact power law, 
p r~j r _/3 , the spherically symmetric and isotropic Jeans 
equations only allow the solutions where the density power 
slope is in the range, 1 < f3 < 3. This result is independent 
of the value of e, and shows that if the central density in- 
deed is a power law, then the density profile cannot be more 
shallow than (3 = 1. 

This constraint weakens slightly for a more general den- 
sity profile where the density power slope, (3{r) , is a function 
of radius. The inner density profile is then allowed to be as 
shallow as described in eq. (12). Also for anisotropic systems 
more shallow profiles are allowed, according to eq. (8) . 
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